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ABSTRACT 

The  influence  of  eroeivity  on  iht  acoustic  stability  of  solid 
fuel  rocket 8  Is  examined.  Allow  ..me  e  ia  made  for  a  frequency- 
dependent  erosion  c  cm  slant.  The  theory  verifies  the  stabilising 
iaflueisce  oi  the  steady  flow  erosivity  for  propellants  having 
positive  erosion  constants,  with  the  contrary  prediction  for 
those  having  negative  erosion  constants.  The  influence  of 
acoustic  erosion  is  examined  and  th  orde"  in  v  oh  it  contributes 
is  shown  to  depend  on  the  symmetry  of  the  flow  field  in  the  rocket, 
'being  quite  different  for  rocks te  with  nozzles  at.  the  end  as  com¬ 
pared  to  rockets  with  nobles  in  the  center  Degenerate  and 
quasi -degenerate  cases  where  erosivity  might  become  a  predomi¬ 
nant  effect  are  discussed. 
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I.  INTRODUCTION 


2 


A  burning  pn  oellant  .  /v  possesses  the  ability  of  amplifying  acoustic 
disturbances  at  its  surface  over  a  relatively  wide  frequency  band.  If  the  net 
acoustic  losses  in  the  rocket  motor  are  sufficiently  small,  as  they  frequently 
are  near  the  resonant  frequencies  of  the  rocket  cavity,  the  motor  acts  as  an 
acoustic  oscillator.  The  high  amplitude  pressure  oscillations  which  then  result 
produce  a  wide  variety  of  secondary  phenomena.  * 

Empirical  studies  have  led  to  a  few  general  observations,  among  which 
is  the  very  interesting  one  that  it  is  the  fluctuation  in  gas  pressure  (in  contrast 
to  gas  velocity)  which  appears  to  L>«_  of  prime  consequence  in  determining  the 
stability.  Corresponding  to  this  observation,  many  theoreticai  approaches  have 
ignored  the  possible  effects  of  erosive  velocity  on  burning  rate  and  have  actually 
achieved  a  considerable  progress  toward  understanding  the  phenomenon.  These 
approaches  are  concerned  chiefly  with  the  acoustic  stability  of  the  burning  pro¬ 
pellant.  The  acoustic  stability  is  determiiutrc!  by  whether  or  not  a  pressure  wave 
(of  some  arbitrary  frequency)  incident  on  the  burning  surface  alters  the  burning 
rate  in  such  a  way  as  to  inhibit  or  reinforce  the  incident  wave  The  boundary 
conditions  at  the  burning  propellant  surface  then  are  expressible  in  the  usual 
admittance  form  and  the  motion  of  the  gas  is  simply  described  in  terms  of  the 
usual  acoustic  modes  of  the  system.  Ultimately,  of  course,  the  acoustic 

^See,  for  example,  ’Analysis  of  Results  of  Combustion  Research  on 
Stolid  Propellants",  E.  W.  Price  Solid  Propellant  Rocket  Research  Confer¬ 
ence,  Princeton  University,  Princeton,  New  Jersey,  Jan,  28-29,  i960, 

APIS  preprint  No.  i CSS -60 
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stability  rests  upon  the  sign  of  the  real  part  of  the  acoustic  admittance,  which 
is  the  indicator  of  amplification  or  attenuation  of  a  sound  wave  at  a  surface. 

A  variety  of  phenomena  now  seem  to  be  accounted  for  in  terms  of  such 

.  2, 3, 4, 5 

analyses.  ’  ’  ’ 

In  spite  of  the  empirical  observation  of  the  prime  importance  of  the 
acoustic  pressure,  and  the  success  of  theories  based  thereon,  one  cannot  help 

2  .  - 

’’Combustion  Instability  Acoustic  Interaction  with  a  Burning  Surface”, 
R.  W.  Hart  and  F,  T.  McClure,  J.  Chem  Phys.  3CP,  1501  (1959).  Also  APL 
report  TG-309  "Effect  of  Solid  Compressibility  on  Combustion  Instability”. 

J.  F.  Bird,  L.  Haar,  R,  W.  Hart,  and  F  T  McClure,  scheduled  for  April 
issue,  J.  Chem.  Phys.  Also  APL  report  TG-335 -1 . 

3 

”On  Acoustic  Resonance  in  Solid  Propellant  Rockets”,  F,  T.  McClure, 
R.  W.  Hart,  and  J.  F.  Bird,  scheduled  for  May  issue,  J„  Appl.  Phys.  A’so 
APL  report  TG-335-2,  "Solid  Propellant  Rocket  Motors  as  Acoustic  Oscillators” 
F.  T.  McClure,  R.  W.  Hart,  and  J.  F.  Bird  Solid  PropellaJtt  .Rocket  Research 
Conference,  Prin^eton  University,  Princeton,  New  Jersey,  Jan,  28-29,  1960, 
ARS  preprint  No.  1049-60  Also  APL  Report  TG-335  3 

4  „ 

"High  Frequency  Combustion  Instability  in  Solid  Propellant  Rockets. 
Part  I”.  Sin -I  Cheng.  Jet  Propulsion  24  27-32  '1954). 

5 

"Unstable  Burning  Phenomenon  in  Double -base  Propellants”.  Theos 
A.  Angelus,  Solid  Propellant  Rocket  Researc  h  Conference,  Princeton  University, 
Princeton,  New  Jersey,  Jan.  28-29,  1960  Also  Allegany  Ballistics  Laboratory 
Report  AEL  Z-9. 
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velocity  might  not  be  important  to  the  stability,  at  least 


under  some  circumstances.  Since  the  steady  state  burning  rate  depends  on 
both  the  pressure  and  the  erosive  velocity  (the  gas  stream  velocity  component 
parallel  to  the  burning  surface)  the  theoretical  basis  of  the  empirical  observa¬ 
tion  and  its  probable  generality  deserves  elucidation. 

Although  the  role  of  acoustic  erosive  velocity  has  not  been  entirely 

g 

ignored,  for  the  most  part  the  analyses  in  which  it  plays  a  part  have  been 
directed  toward  the  question  of  inherent,  rather  than  acoustic  instability. 
Inherent  stability-  is  concerned  with  whether  or  not  an  arbitrary  perturbation 
on  the  burning  rate  decays,  rather  than  grows,  exponentially  with  time.  A 
’’propellant”  which  is  inherently  unstable  would  not  be  expected  to  possess  a 
’’steady  state”  of  combustion  under  any  ordinary  circumstances.  On  the  other 
hand,  we  recognize  that  a  propellant  for  which  transients  are  damped,  i.e.  . 
which  possesses  a  steady  state,  may  or  may  not  act  as  an  acoustic  amplifier 
for  incident  sound  waves. 

The  bearing  of  erosive  velocity  on  acoustic  stability  was  discussed 

ry 

briefly  by  Cheng  who  concluded  that  "the  dependence  of  the  burning  rate  on 
the  drifting  velocity  alone  is  not  likely  to  excite  instability.  .  . "  Since  the  scope 
of  his  analysis  was  rather  limited,  this  conclusion  cannot  be  accepted  as 


definitive. 


uSee,  for  example.  L.  Green  and  W  Nachbar.  "Analysis  of  a 
Simplified  Model  of  Solid  Propellant  Resonant  Burning".  J.  Aero.  Space 
Sc iencto  26.  518  (1959). 

(Sin-I  Cheng.  "High- Frequency  Combustion  Instability  in  bolid  Pro¬ 
pellant  Rockets."  Part  J,  (particularly  Appendix  II).  Jet  Propulsion  24, 
102-109  (1954). 
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ine  suuject  which  we  shall  consider,  tnen,  is  tnai  oi  acoustic  stauuuy 


in  the  presence  of  a  burning  rate  dependent  on  both  pressure  and  erosive 
velocity.  We  intend  to  ask  for  the  conditions  under  which  such  a  burning  sur¬ 
face  amplifies  or  attenuates  an  incident  sound  wave.  The  question  of  stability 
against  finite  disturbances  is  auite  another  problem,  which  will  not  be  treated 
here. 


II.  ANALYTICAL  FORMULATION 


In  order  to  suppress  unnecessarily  complicating  distractions,  we  shai, 
simplify  .he  geometry  to  that  of  a  circular  cylindrical  cavity  containing  a 
cylindrical  grain  (outer  radius  b,  inner  radius  a).  We  shall  neglect  the  small 
doppler  effect  produced  by  the  rather  slow  mean  flow  of  burned  gas  and  regard 
the  burned  gases  as  a  homogeneous  medium  insofar  as  acoustic  motion  is  con¬ 
cerned.  For  our  present  study,  which  is  concerned  merely  with  the  relative 
importance  of  pressure  and  erosive  velocity,  we  will  select  a  method  by  which 
we  may  also  neglect  acoustic  losses  in  the  gas  and  in  the  nozzle. 

2  I  The  End  Face  Boundary  Conditions 


We  shall  treat  only  the  case  of  rigid,  lossless  ends  so  that  we  have 


;Vz.=e =  0 


and  (u  )  = 0  . 

z  z  --  c, 


(1 


i  i  a  i 


where  z  -  0  and  z  =  L  are  the  coordinates  of  the  end  faces  in  a  cylindrical 
svstem  whose  axis  coincides  wun  that  of  our  cylindrical  cavity,  and  u_  is  the 
acoustic  velocity.  These  boundary  conditions  correspond  to  those  encountered 


in 


some  researen  type  rocket  motors  with  sip.  nozzles,  and  may  be  thought  to 


approximate  end  nozzle  rockets  without  end  cavities  and  with  high  port  tt 


throat  rat  it's. 
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The  formulation  of  th^  erosive  boundary  condition  at  the  propellant  sur¬ 
face  is,  of  course,  the  rrux  of  the  problem,  and  at  the  very  beginning  we  are 
somewhat  nonplussed  by  the  absence  of  a  time -dependent  theory  of  the  erosion 
of  the  burning  boundary  Laver.  We  may  elect  to  proceed  in  either  of  two  direc 
tions.  On  the  one  hand,  we  might  elect  to  stop  and  attack  the  development  of 
such  a  theory,  or  on  the  other  hand,  we  might  attempt  to  formulate  the  acoustic 
analysis  in  sufficiently  general  'erms  that  the  resuits  of  a  future  theory  could  be 
inserted  into  the  analytical  structure.  These  two  paths  are  rather  distinct.  The 
former  course  would  apparently  involve  a  conglomeration  of  chemistry  and 
aerodynamics  whereas  the  latter  is  primarily  an  application  of  acoustical 
principles.  Rather  than  concern  ourselves  with  the  predictions  of  a  particular 
burning  theory,  we  elect  to  oursue  the  latter  course  and  focus  our  att°ntion  on 
the  development  of  a  rather  general  acoustic  structure.  By  following  this  path, 
we  incidentally  may  hope  to  be  able  to  form  an  estimate  as  tc  the  importance  of 


a  time -dependent  erosive  burning  theory. 

The  erosive  boundary  condition  may  be  obtained  from  the  general 
relationship  expressing  the  flux  density  of  hot  gas  generation,  in.  Restricting 
our  attention  to  the  periodic  case  -i.  e.  .  the  transients  have  already  decayed), 
we  may  express  n.  as  a  function  of  pressure,  the  erosive  velocity,  and  their 

m 

time  derivatives  Thus. 


m  ic  nn  P.o  .  ;d  ’  ... 
a  ‘at 
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Lnte<  cals  wah  respect  to 


me  would  also  be  required  in  the 


rails  lent 


case. 
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surface  (r-a),  vo  is  the  erosive  velocity  (i.e.  ,  gas  velocity  parallel  to  burning 
surface)  and  where  the  subscript  t’s  indicate  partial  differentiation  with  res¬ 
pect  to  time.  The  erosive  burning  depends  on  the  absolute  magnitude  of  v 
rather  than  on  v  itself,  because  the  effect  of  erosion  is  presumed  to  be 
independent  of  the  direction  of  the  erosive  velocity. 

2.3  Perturbation  Form  of  the  Erosive  Boundary  Condition 

i  | 

Assuming  the  burning  rate  to  be  an  analytic  function  of  j  |  ,  its 
time  derivatives,  and  p  and  its  time  derivatives,  we  may  expand  about  the 
point  (P;  0.0,0.  .  .  .  ;0. 0,0,0.  . .  ),  i.e.  ,  about  the  non -os<  dilatory  stale  with  no 


erosion.  Thus,  we  have 

m  =  [m(P  :0. 0, . . . )  *  <  O- 

L  /  CPa 
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+  V  higher  order  non -linear  temr  »  !  i 

L  _  j  J 

where,  to  avoid  urwieldy  notation,  we  have  displayed  only  one  of  tne  multitude 
of  subscripts  "zero”  attached  to  pach  of  the  partial  derivatives.  Insofar  as 
the  onset  of  instability  is  concerned,  we  may  neglect  higher  order  terms 
provided  the  steady  state  erosive  burring  law  •  ;  adequately  approximat'd  by  a 


br»m  r  tun  'tuvn  t  v 

until  I  iU.‘U.  iiUu  v  i.  v 


The  first  term  on  the  rhs  Eq.  1c;  .s  mere  :v  the  dc  burning  rate 


the  absence  of  erosion,  which  we  designate  by  nip  .  The  first  of  tta 
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bracketed  sums  expresses  the  contribution  to  gas  flux  of  a  fluctuating  pressure 
with  no  erosion.  In  the  notation  of  ref.  2.  this  bracket  will  be  denoted  by 
(mQ  The  second  bracket  on  the  rhs  of  Eq,  (lc)  describes  the  effect  of  an 

erosive  velocity,  with  no  pressure  fluctuation.  (We  may  recognize  the  term 

-  as  uif.k  ,  where  k  is  the  usual  (dc.)  erosion  constant.)  The  impor- 

8  v  ift 
i  e  0 

tant  point  here  is  that  the  pressure  and  erosive  contributions  to  burning  rate 
are  now  separated,  with  the  erosive  part  expressed  as  a  linear  function  of  vg 
and  its  time  derivatives.  The  pressure  part  has,  of  course,  been  studied 
separately,  elsewhere.  * 


Finally,  we  complete  the  perturbation  formulation  by  recalling  that  m  may 
be  expressed  as 


-m 


v  ) 

s  r-a 


(If) 


where  p  is  the  hot  gas  density,  v  is  the  hot  gas  radial  velocity,  and  v  the 

regression  spe-.d  of  the  solid  surface,  with  radial  velocities  measured  positive 

in  the  outward  radial  direction.  Wc  have 

v  =  v  +  Y  p 
s  s  s  a 

wh‘jre  the  bar  denotes  the  time  average,  Y  -  specific  acoustic  admittance  of 
the  solid  surface,  which  is  zero  it  the  solid  is  rigid,  and  -  the  acoustic 
pressure  ai  the  propellant  surface,  where  r  ---  a.  We  use 


P  P  1  f 


y  p 


dg) 


where  v  is  the  hot  gas  specific  heat  ratio  and  P  p  —  is  (he  mean 
‘  -  > 
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chamber  pressure,  with  c  the  sound  velocity,  p  the  mean  density  of  product 
gases.  Equations  (If)  and  (lg),  with  the  definitions  of  mg  and  mg/ij  given 
above,  when  substituted  into  Eq.  (lc),  yield 


For  the  present  case,  one  further  simplification  of  the  boundary  equation 
also  applies.  Noting  that  the  magnitude  of  the  erosive  velocity  may  be  expressed 
in  terms  of  its  (cylindrical  coordinate)  components  as 


ve|  *V  l  ve  +  uz>  + V 


we  expand  the  radical  to  see  that,  for  the  arbitrarily  small  acoustic  amplitudes 
relevant  to  stability  theory, 


e 


A / 

\  S\J 


V  +  u 
e  z 


=  ve'|  +  uz  Signum  (ve), 


(li) 


where  Signum  (v  )  =  +  )  or  -  l  according  to  whether  vo  is  positive  or  negative. 


In  the  following,  then,  we  will  use  the  linear  boundary  condition  obtained  by 
substituting  Eq.  (li)  into  Eq.  (Ih),  i.e., 


/C 


u  \  pcY  t 

j  -  M  (p.  -  — )  + - —  -  -  M  k 

c  /  avl  y  y  a 

y  r=a 


—  I  €  i 

V  j - + - 

c  y 

pc 


Signum  (v  ) 


•  («jt 
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3  m 


(2) 


10 


in  which  v  is  still  to  be  regarded  as  an  arbitrary  funcf  n  of  z  . 

2.  4  Fourier  Transformation 

It  is  most  expedient  to  take  advantage  of  the  simplifications  which  result 
when  we  consider  only  harmonic  solutions.  Accordingly,  we  may  define  the 
time -independent  tilda'd  functions  by  the  relationships 

p  =  (Real  of)  'p 

u  =  (Real  of)jue*w*  ,  and 

=  (Real  of)  ,  (2a) 

The  boundary  equation  (Eq.  (2)  )  then  takes  the  form 

M  k  fv  1 

e  =  - - +  M  K  u  Si  gnu  m  (v  ) 

y  a  z  6  e 

;2b) 


where  K,  because  of  this  definition,  may  readily  be  obtained  from  a  time 
dependent  theory  of  erosion  if  it  should  become  available.  We  may  note,  in 
particular,  that  0)—>  k,  and  refer  to  K  as  the  ”ac  erosion  constant”. 

The  end  face  boundary  conditions  (Eq.  (la)  )  are  easily  transformed  to  read 

= 0  (2c) 

Z’z=0  z=L 


and,  of  course,  the  wave  equation  to  be  solved  for  these  boundary  conditions 
transforms  to  the  Helmholz  equation,  'p  ^(^~—\  'p  =  0  . 

Vc  / 

III.  ANALYTICAL  SOLUTION  OF  THE  BOUNDARY  VALUE  PROBLEM 
The  analytical  problem  has  now  been  defined,  ff  the  solution  were  obtain¬ 
able  by  several  of  the  mere  obvious  and  routine  approaches,  we  would  merely 
display  the  result,  and  assert  that  it  easily  could  be  verified  by  substitution 
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into  the  Helmholz  and  the  boundary  equations,  A  rather  special  approach  seems 
necessary,  however,  so  that  it  seems  desirable  briefly  to  sketch  the  method  of 
solution.  The  method  begins  in  the  standard  way.  We  note  that  the  solution 
will  be  expressible  as  a  sum  of  the  elementary  (separable)  solutions  of  the 
Helmholz  equation  expressed  in  cylindrical  coordinates: 

CO  00 


PI.J  Jj(alr>  cos(% 


cos  (j0) 


(8a) 


00  00 


U 
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-c 
i  yw 


in  z 


pf  ja|  Jj(a|r)  cos  (  — -  !  cos(j0) 


(3b) 


ao  oo 
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r~> 

u  * - 

z 


ivz 


iyo;L 

j  =0  1  =  0 


pf  .  I  J.(d1r)  sin  /  — cos  (j0) 


(3c) 


The  appearance  of  the  individual  sine  and  cosine  terms  indicates  that  we  have 
forseen  the  axial  bound  a*  y  conditions  expressed  by  Eq.  (2c). 

The  usual  approach  would  ue  to  substitute  the  field  expression  (Eqs.  (3a) , 
(3b),  (3c))into  the  boundary  equation  (Eq.  (2b)  ).  In  the  absence  of  erosion 
(k  -  K  =  0),  this  substitution  quickly  yields  the  usual  result  that  the  charac¬ 
teristic  frequencies  would  be  determined  by  the  values  of  u>  for  which 
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)  y  a  <£ 


■(a  a)J,(ot  a)  +  - 

v  n  '  j'  n  c 


pc  Y 


S 


~yK 


Ml. (a  a) 
a  jv  n  ' 


We  note  that  the  system  would  be  neutrally  stable  (real  o>)  only  if 


Real  of 


=  -45-  Heal  of  Y  , 
yU  s 


a 


and  that  the  system  would  be  damped  only  if 

i'\ 

m 


Real 


.  Hcal  of  Y  . 

y  7  y  Ma  s 
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=  0 


(4a) 


(4b) 


In  the  presence  of  erosion,  however,  the  analytical  difficulties  begin 
immediately.  Commencing  as  before,  we  note  that  in  order  to  determine  the 
characteristic  frequencies  (for  which  a  solution  for  the  p^'s  is  at  least  possible), 
we  would  have  to  discover  the  roots  of  an  exceedingly  intractable  infinite  determi¬ 


nant. 


At  this  point,  a  bit  of  reflection  is  enough  to  convince  us  that  we  did  not 
really  want  specifically  to  know  about  the  effect  of  erosion  on  the  characteristic 
frequencies,  anyway.  The  information  which  we  really  seek  is  the  relative 
importance  of  erosion  compared  to  pressure  in  determining  stability. 

Let  us,  then,  attempt  a  less  standard,  but  more  direct  approach.  Suppose 
we  recognize  that  we  would  have  a  rather  satisfactory  answer  to  our  question  if 
we  knew  for  a  given  erosion  constant  how  much  the  admittance  (as  measured, 
say,  by^/cf  of  the  burning  surface  would  have  to  be  altered  in  order  to  restore 
the  characteristic  frequency  to  its  value  for  zero  erosion.  The  important 
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advantage  in  formulating  Jiis  particular  question  is  that  the  characteristic 
frequencies  of  the  system  then  correspond  exactly  to  those  for  zero  erosion, 
which  are  very  easy  to  determine! 

In  order  to  discuss  the  contribution  of  erosive  velocity  to  stability,  we 
must  first  digress  for  a  moment  in  order  to  display  the  nor  lenclature  for 
labeling  the  acoustic  modes.  To  specify  a  particular  mc.ie,  ve  must  give  angular, 
radial  and  axial  mode  numbers.  Thus,  a  mode  is  represented  in  the  usual  way 
by  a  three -element  array,  the  first  number  representing  the  number  of  axial 
half  waves,  the  second  the  number  of  radial  interior  (velocity^  .odes,  and  the 
third  the  number  of  angular  nodes  In  order  to  avoid  unnecessary  repetition  of 
subscripts,  we  shall  designate  the  modes  by =  N, R,6. 

In  order  to  carry  through  this  approach,  then,  we  shall  write  k  -  /3K  , 

* 


and  develop  the  solution  as  a  power  series  in  K.  We  write,  for  the^Nth 

mode, 

7 


j  =  1 


Kc  a  M 


(where  the  expansion  parameter 


,  rather  than  merely  K. 


itself,  is  used  for  later  convenience).  Thus,  if  we  can  determine  the 


* 

The  value  of  /I  would  have  to  be  determined  from  a  time-de'iendent 
erosion  theory  and  in  general  would  be  a  function  of  propellant  composition, 
operating  conditions  and  frequency.  If  there  were  frequency  intervals  for  which 
k  was  much  greater  than  K,  the  appropriate  expansion  for  our  problem  would 
of  course  be  in  powers  of  k,  with  K  =  fi'k 


I  l 
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coefficients  D,  VT  we  shall  have  found  the  comparison  between  the  effects  of 
j,  N 

/yw 

erosive  velocity  and  pressure.  Restricting  our  attention  to  the  conditions  for 
neutral  stability,  for  which  cu's  are  real  (and  given  by  Eq.  (4a)  ),  one  can 
now  readily  verify  that  the  solution  to  our  boundary  value  problem  is  expressed 
by  Eq.  (3a)  with  the  given  for  the^th  mode  of  the  system  by 


f  aKcM  \n  {n\ 

=  6  >  ( _ _  |  p y  N  ,  (5a) 

■J8.  i'e  A  V  L  /  1  ~ 

n  =0  7 


t3 


where  the  _T  are  found  from  the  recursion  formula 

*  l  ,N 


port  mach  number. 


port  mach  number. 


for  If ■  N 
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This  completes  the  solution  for  the  non-degenerate  case.  Calculation  is 

not  difficult  in  spite  of  the  rather  unwieldy  expressions  because  of  the  typi - 
,  f  aKcMa  \ 

cally  small  value  of  |  - - — - \  .  It  will  usually  be  entirely  adequate  to 


\  ' 

(2)  ' 

terminate  with  the  pv  '’s  In  a  degenerate  case,  where  more  than  one  mode  of 
motion  corresponds  to  a  single  frequency,  the  denominator  of  Eq.  (5b)  vanishes 
for  some  If  N,  and  a  slightly  more  general  soluuon  must  then  be  obtained. 


IV.  RESULTS  OBTAINED  FROM  THE  ANALYTICAL  SOLUTION  OF  BOUNDARY 

VALUE  PROBLEM  (NON- DEGENERATE  CASE/ 


The  question  to  lie  considered  in  this  section  is  the  effect  of  an  erosive 

velocity  on  stability.  As  has  been  discussed,  we  have  decided  to  measure  this 

effect  by  determining  the  increment  in  the  response  function  (p/€j  which  would 

be  necessary  to  re-establish  the  neutral  stability  condition,  when  erosion  is 

suddenly  ’’turned  on”.  The  reference  scale  is  provided  by  the  observation  that 

the  real  part  of  p/ e~  can  ordinarily  be  expected  to  lie  between  the  dc  pressure 

2 

exponent  (<C1)  and  a  value  perhaps  several  tenths  greater  than  unity.  Of  course, 
we  must  keep  in  mind  that  the  quest iun  of  instability  or  stability  is  determined  by 


whether  the  real  part  of  p  e  does  or  does  not  exceed  a  critical  value,  and 


therefore  the  erosive  term,  though  small,  would  n 


vertheless  become  a  erpical 


factor  in  determining  stability 


for  eases  such  that 


<)q  were  close  to  the 


critical  value. 

Examining  the  magnitude  < 

we  suspect  that,  largely  because 
hot  gases  leaving  the  burning  zo: 


f  the  expansion  parameter 
oft  he  s  ;na  I  i  mac  h  nu  mbe  r 


i 


a  M  Kc 
a 

L 


/ 


/ 

D 


u  s 


ually  will  be  quite  satisfactory  u 


17 


consider  only  the  lowest  order  non-vanishing  terms  even  for  highly  erosive 
propellants. 

Let  us  consider  the  linear  term  first.  Substituting  Eq.  (5d)  into  Eq. 
(5c),  and  the  result  into  Eq.  (4c),  we  find  rather  quickly  that  (we  also  use  the 
identity  0K=k), 


=  — —  M  bxrkT  +  c.  VT  +  higher  order  terms.  (6) 

y  p  NN  ya'L  N,  N 

0 

In  discussing  Eq.  (6)  two  cases  need  to  be  distinguished,  namely. 

the  ordinary  rocket  with  the  nozzle  at  the  end  and  special  rockets  with  a  side 

nozzle  at  the  center.  In  the  first  case  v^  is  a  positive  quantity  whereas  in  the 

second  case  it  is  positive  in  one  end  of  the  rocket  and  negative  in  t..e  other. 

For  the  usual  case  of  positive  v  the  coefficient  c.T  XT  =  0.  Thus,  to 

e  N.  N 

first  order,  the  contribution  of  the  acoustic  erosion  is  zero.  For  the  typical 

-  1 

cylindrical  rocket  v  —  M  cz,  L  we  obtain  bvrXT  .  Tnus.  we  have  already 

e  p  NN  2 

determined  in  this  first  order  term  the  predominant  effort  of  mean  erosive 
velocity  on  the  onset  of  instability.  It  is  clear  that,  with  kc  typically  less  than 
about  4  even  for  very  erosive  propellants,  and  M  typically  not  greater  than 


per  naps  1.  10.  the  effect  of  dc  erosion  is  to  be  "compensated  for"  by  altering 
a  e  by ^  1  5.  Thus,  we  should  expect  to  find  enhanced  stability  for  pro¬ 
pellants  having  a  positive  dc  erosion  constant .  For  unusually  high  erosive 
velocity,  and  mghiy  erosive  propellants,  this  stabilizing  effect  could  be  very 


substantial.  Such  an  effect  has  indeed  be^n  observed.  On  the  other  hand,  we 


would  expect  t  he  opposite  e fleet  for 
greater  tendency  toward  instability . 


negative  e.  os  ion  constants,  namely. 
In  under  to  discuss  the  effect  of  the 
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acoustic  velocity  for  the  end  nozzle  case  we  must  go  to  second  order.  This 
we  will  do  after  discussion  of  the  first  order  treatment  for  the  center  nozzle 
case. 


For  the  center  nozzle  rocket  the  situation  is  quite  different.  Here 

L 


— 

v  - 
e 


2M  cz 
P 


z  < 


V 


sU 
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2M  c(L-z) 
P 
L 


This  steady  state  velocity  distribution  leads  to 

2 

c  VT  .  =  -rr -  for  N  odd 

N,  N  Ntt 

c  vt  VT  -  0  for  N  even; 

N ,  N 


while 

h  1  2 

NN  '  2  '  ~2~2 
v  N 

1 

=  2 


for  N  odd 


for  N  even 


Thus  we  see  that  for  N  even  (e.g.  ,  the  even  axial  modes!  the  cc-nter  nozzle 
rocket  gives  the  same  results  as  the  end  nozzle  rocket.  However,  for  N  odd 
{e.g,  ,  the  odd  axial  mode  '  the  situation  is  strikingly  different .  Here  the 
contribution  of  the  acoustic  erosion  could  be  extremely  important.  In  oar- 


ticular.  for  the  first  axial  mod*,  we  find 


\ 


k : ' 


o 

"K 


2  i  Kc 


t  /. 


and  we  note  that  rhe  ordinarily  stab:  dz.oc  i for  k  0  effect  of  the  steady 


IS 


erosion  is  reduced  by  40%.  More  important,  the  imaginary  part  of  K  need 
not  be  large  to  make  an  extremely  significant  contribution.  For  example, 
if  imag  (Kc)-  -  i/2,  the  acoustic  term  would  be  1  />  ,  which  is  the  usual 
criterion  on  (jl/f€)  for  an  amplifying  surface 

Provided  k  is  positive,  the  damping  due  to  the  b^  ^  term  is  greater  for 
even  modes  than  for  odd  modes.  The  first  order  acoustic  term  contributes  to 
the  instability  if  (iK)<£  0,  and  to  the  damping  if  (iK)>0.  Thus,  for  low  frequency 

(and  k>0)  one  would  expect  this  term  to  contribute  to  instability  ior  tne  odd 

* 

modes.  One  cannot  help  but  w^  ,der  whether  this  effect  does  not  contr  bute  to 
the  experimental  obser  ation  that  the  predominantly  unstable  axial  modes  m  the 
center  nozzle  of  Price  et  al  are  the  odd  modes.  However,  at  much  higher 
frequencies  the  odd  modes  might  be,  on  the  contrary,  stabilized. 

Returning  to  the  discussion  of  the  end  nozzle  case,  we  shall  investigate  the 
contribution  of  acoustic  erosion  by  displaying  the  second  order  terms.  We  find 


v,  l 
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amplitude  because  the 
itself,  would  then  gene 


do  not  pertain  to  the  established  oscillation  at  high 
mn - . meant les ,  particularly  those  due  to  erosion 
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and 


b 


£,£ 


1 
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Substitution  of  these  expressions  into  the  general  solution  finally  yields,  for 
this  case,  n 

kcM 


/-V- 
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/V-1 


ii_  _  _  p 
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While  Eq.  (7)  may  appea.  simpler  than  the  exact  solution,  this  appearance 
is  SGmewhat  deceptive  and  one  must  be  careful  about  generalizing  from  it  as 
to  importance  of  the  erosive  effects.  These  obviously  depend  upon  the  nature 
of  the  propellant  (particularly  with  respect  to  the  values  of  its  erosion  constants 
K  and  k)  and  upon  the  geometry  and  the  mode  under  consideration.  If,  for 
example,  the  denominator  in  the  sum  should  fortuitously  become  very  small 
for  some  particular  £  and  J^the  result  would  be  a  large  effect  on  the  stability, 
either  to  increase  or  decrease  it  depending  on  the  sign  of  the  real  part  of  the 
term  involved.  The  extreme  case  of  one  of  these  denominators  actually  being 


zero  corresponds  to  the  degenerate  case  discussed  in  the  section  V. 


4 
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In  many  typical  rocket  geometries,  however,  the  denominators  will  be 

more  or  less  uniform  in  v'alue  over  the  frequency  range  of  interest  and  in  that 

event  one  can  expect  less  drastic  results.  In  the  case  of  N  =  N,0,Q,  or 

a-  n  L 

essentially  pure  axial  modes,  we  nc'ice  that  — -  =  n  and  as  a  conse- 

Nc 

quence,  for  the  usual  motors  and  propellants,  the  second  term  in  each  square 

bracket  (k  terms)  will  make  a  negligible  contribution  if  the  axial  number  N  is 

not  very  large.  The  contribution  of  the  first  terms  (K  terms),  however,  needs 

examination  since  they  may  have  larger  coefficients.  To  gain  some  idea  of 

the  size  of  the  effect  let  us  consider  the  special  case  with 

(m/cJ0  -  (l/y)  -  {pcY  g/yMa)J  =  0  which  corresponds  to  the  unperturbed  soiu 

tion  (see  Eq.  (4a))  having  a  pressure  antinode  at  the  surface.  Assuming  the 

rocket  to  be  a  cylinder  with  v  =  M  cz/'L  we  obtain 

e  p 


M 


**  \  k^M 
Jf_  ]  d  P 


e 


ST  +  (1 4  V  o> 


-1 


2iM  (Kc)' 

3. 


Y* 


2  r. 


v 


1-  (-D 


i/+N 


'.(r/rw) 


,  2  XT2,2  f.  2  /x  2  _  N?ra  /T  TZZ. 
{v  -N  )  a/  1  -  v  J N  J jV — x7”V  1  "  v  /N 


(C  a) 


N  7TR 

When  — ^ —  is  small,  the  sum  converges  rather  rapidly  because  of  the  high 
power  in  the  denominator,  and  clearly,  the  largest  contributions  will  ordi¬ 
narily  come  from  the  terms  v  -  N  +  1,  except  for  the  case  N  =  1  when  only 


the  v  =  N  +  1  terms  need  be  retained.  In  particular,  for  the  lowest  axial  mode; 

LM_  M 

we  find  (with  N  =  1,  ©  0,  and  v  -  2  and  - 
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For  this  particular  case  the  acoustic  contribution  is  proportional  to  the 
imaginary  part  of  K//3  (i.e. ,  (K^/k)  )  so  that  there  is  no  contribution  to 
stability  unless  K  is  complex.  The  importance  of  the  result  depends,  of  course, 
on  the  size  of  this  component  of  K/0.  If  it  should  be  of  the  sa^e  order  as  k 
itse^,  i.e. ,  kc  say  about  4,  the  acoustic  contribution  would  be  approximately 
one -half  the  dc  contribution.  Thus  for  positive  erosion  the  net  effect  would  be 
to  stabilize  the  rocket. 

Investigation  of  Eq.  (7)  for  radial  and  tangential  modes  shows  that  the 
acoustic  term  is  quite  small  for  typical  geometries.  Thus,  we  see,  in  con¬ 
trast  to  the  center  nozzle  rocket,  there  is  a  wide  range  of  conditions  over  which 
we  expect  erosivity  to  contribute  little  to  the  question  of  stability  or  instability 
For  rockets  in  which  axial  modes  of  very  high  N  are  not  heavily  damped,  this 
generalization  may  not  be  valid.  Moreover,  it  may  fail  dramatically  for 
conditions  approaching  the  degenerate  case  discussed  in  the  next  Section. 

V.  THE  DEGENERATE  CASE 
5. 1  Degeneracy  and  Quasi -Degeneracy 

If  we  were  to  consider  maximizing  the  acoustic  erosive  effect,  we  would 
refei  to  Eq.  (7),  and  note  that  if  the  denominator  were  to  become  very  small 
for  some  term,  say  for  l  -  M,  then  the  influence  of  erosion  might  indeed 
become  overwhelming.  Of  course,  the  real  part  of  the  denominator  will  vanish 
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whenever  the  frequencies  of  two  characteristic  mod^s  of  the  system  coincide. 
Since  the  internal  radius  of  the  usual  solid  propellant  rocket  chamber  increases 
with  time  as  the  propellant  burns,  the  frequencies  of  the  radial  and  tangential 
modes  decrease  with  time  while  the  pure  axial  mode  frequencies  are 

independent  of  chamber  radius.  This  state  of  affairs  is  indicated  in  Figs.  1 

5  * ** 

and  2,  for  a  hypothetical  motor  with  L  =  5rcm,  c  =  10  cm/sec  .  In  Fig.  1 
we  have  plotted  the  radial  and  axial  frequencies,  and  in  Fig.  2,  the  first 
tangential  frequencies,  corresponding  to  the  radial  velocity  node  boundary 
condition  on  the  burning  surface. 

In  general,  at  each  intersection  point  in  Figs.  1  and  2,  the  real  part  of 
the  denominator  vanishes  in  one  of  the  p^'s  occurring  in  the  sum  representa¬ 
tion  of  the  solution.  But  the  system  is  not  truly  degenerate  at  these  intersec¬ 
tions  unless  the  imaginary  part  of  the  denominator  also  vanishes  there,  i.  e. , 


* 

The  figure  is  easily  scaled  to  represent  chambers  of  other  lengths  and 
radii.  If  both  a  and  L  are  multiplied  by  a  common  factor,  the  frequency 
scale  is  to  be  divided  by  that  factor. 

** 

Although  these  "quasi-modes”  are  not  the  true  modes  of  the  cavity 
because  the  solid  propellant  itself  also  indulges  in  the  acoustic  motion,  it  has 
been  shown  elsewhere  that  the  frequencies  of  the  modes  of  the  composite 
system  whi,ch  have  large  gas  pressure  amplitude  -  and  therefore  large 
erosive  acoustic  velocity,  cf.  (Eq.  (3c)  )  at  the  propellant  surface  lie  close 
to  these  gas  quasi -mode  frequencies.  3 
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Fig.  1  FREQUENCY  VERSUS  GEOMETRY  FOR  THE  AXIALLY  SYMMETRIC 
MODES  OF  A  GAS  CONTAINED  IN  A  UIG ID-WALLED 
CYLINDER  (GAS  QUASI-MODES) 

The  axial ,  radial,  and  a.-^ulur  modi’  numbers  are  shown  for 

each  curve,  m  that  order. 
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Fig .  2  FREQUENCE  VERSUS  GEOMETRY  FOR  THE  FIRST  AZIMUTHAL 
MODES  OF  A  GAS  CONTAINED  IN  A  R I G I D-WALEED 
CYLINDER  (GAS  QUASI -MODES) 

The  axial,  radial,  and  angular  mode  numbers  are  shown  for 

each  rurvi ,  in  that  order . 
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unless  the  acoustic  losses  '’accidentally"  balance  the  gains  for  that  particular 
geometry.  However,  we  see  that  during  the  course  of  burning,  many  points  of 
"quasi -degeneracy"  may  appear,  where  the  denominator  of  a  term  in  the  sum  is 
especially  small  because  its  real  part  vanishes.  But  we  shall  expect  a  true 
degeneracy, characterized  by  the  vanishing  of  both  real  and  imaginary  parts,  to 
be  a  rare  occurrence. 

Unfortunately,  at  this  time  it  seems  impossible  to  calculate  the  expected 
effect  on  stability  at  the  points  of  quasi-degeneracy.  The  impasse  arises  because 
we  need  to  know  the  imaginary  part  of  the  denominator  which  will,  in  general, 
be  a  function  of  the  losses  in  the  normal  mode  being  considered.  The  detailed 
computation  of  these  losses  poses  formidable  difficulties  and  no  experimental 
information  bearing  directly  on  this  point  is  available. 

We  may, however,  expect  to  gain  insight  as  to  whether  or  not  erosion  can 
have  a  large  influence  at  configurations  corresponding  to  such  quasi -degenera¬ 
cies  by  investigating  the  limiting  case  of  true  degeneracy. 

5.  2  The  Degenerate  Case 

Let  us  now  suppose  that  two  normal  modes  of  the  system  are  neutrally 
stable.  In  other  words,  let  the  inner  radius  a  approach  a  value  such  that  the 
denominator  of  Eq.  (7)  vanishes  for  some  l  .  say  f  =  M.  Then,  we  must 
re-solve  the  boundary  value  problem  to  develop  a  solution  for  this  degenerate 
case. 

The  analvtieal  method  remains  relatively  unchanged;  we  again  expand  in 
powers  of  the  expansion  parameter,  and  equate  the  coefficients  of  equal 
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powers  in  the  boundary  equation.  Tins  time,  of  course,  the  unperturbed  motion 
must  include  two.  rather  than  one,  zero  order  terms,  in  order  that  the  new 


sum  to  be  obtained  will  not  contain  any  terms  with  zero  denominators.  We  write 


P 


p  =  PN,N  jJe(aN,Nr)  C0S(^TJ  +  J0<aM,_Nr) 


AJe<aN,Ma)\  Amz'v? 

C08(-x-)tC0S 


f 


♦ 


p<Jj(a|.Nr)  C0S 


(f  M,  N  when  j  =  0) 


cos(jfl) 


(9a) 


where  A  is  to  be  determined  in  such  a  way  that  the  solution  approaches  the 

bracketed  term  as  the  erosion  vanishes  (K-*  0).  The  velocities  and^  are 

r  z 

obtained  in  the  usual  way  (Eqs.  (3b),  (3c)  )  ,  and  substitution  into  the  boundary 
condition  equation  at  r  =  a  finally  yields 


A  = 


0M. 


M,  N 


IW'L  /K  K  \  ! 

(bN.  N  "  bM.  M  i 


:  -f 


+  i 


v 


W 


(9b) 


where 


=  rrNc 


N.  M 


iwL 

c 


N.  M 


(9c) 


The  upper  and  lower  signs  correspond  to  the  two  possible  choices  for  A  in 
Eq.  (9a),  and  define  the  two  orthogonal  combinations  of  the  degenerate  modes, 
which  may  be  obtained  by  substituting  these  A's  into  Eq.  (9a).  The  increment 
in  the  response  function  necessary  to  re-establish  the  neutral  stability  condition 
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when  erosion  is  "turned  on"  is  finally  found  to  he 


kcM  I 

-2T  (bN,  N  >bM.  M >7  (Kc) 


2/c 


u;L  y 


kcM 


l2 

(bN,N  ~bM,M)  l 


Noting  that,  at  least  fo^  typically  small  M^,  (-f^  M  f^  N)  is  positive,  we  see 

that  the  mode  corresponding  to  the  choice  of  the  lower  sign  is  made  more 

stable  by  erosion,  while  the  mode  corresponding  to  the  choice  of  the  upper 

sign  is  made  less  stable.  Each  mode  will  look  more  or  less  like  an  N  mode  or 

Na) 

an  M  mode  according  to  whether  or  not  — - — —  is  substantially  less 

J0(Q?M,Na) 

than,  or  greater  than,  unity. 

An  appreciation  of  the  order  of  magnitude  involved  may  be  obtained  by 

evaluating  Eq.  (9d)  for  a  few  cases.  If  we  consider  the  conventional  end  nozzle 

rocket,  the  b^T  ^  may  be  obtained  from  Eq.  (6b),  and  the  c^  ^  from  Eq.  (6a). 

One  quickly  finds  that  the  increment  in  4/e  arising  from  acoustic  erosion  will 

be  sometimes  large  and  sometimes  small  (relative  to  1/v),  with  the  larger 

contributions  occurring  for  N  close  to  M.  If  we  continue  to  take  Kc  =  4,  L  =  5rr  cm, 
5 

c  -  10  cm,  sec ,  and  >■  =  1.2.  we  find  for  example,  that  for  the  5,2,1  and  6.1,1 
intersections,  the  radical  in  Eq.  (9d)  yields  the  value  1,8.  Thus  such  true 
degeneracies,  if  they  were  to  occur,  could  enable  erosion  to  exert  a  profound 
effect  on  stability.  The  importance  of  .having  M  and  N  nearly  equal  can  be 
shown  by  reducing  the  value  of  M  from  5  to  1  in  the  above  example,  with  the 
consequent  diminishing  of  the  contribution  from  1.  8  to  0.  1. 
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In  the  neighborhood  of  degeneracy  we  see  that  we  might  expect  acoustic 
erosion  sometimes  to  have  a  large  effect  on  the  question  of  stability  or 
instability.  The  question  naturally  arises  why  this  effect  has  not  made  its 
presence  obvious  in  experimental  studies.  Of  °ourse,  one  should  note  that 
unless  the  results  were  relatively  catastrophic,  it  would  be  quite  easy  to  miss 
the  connection  between  degenerate  modes  unless  one  were  specifically  looking 
for  it.  However,  there  are  probably  mure  important  reasons  why  the  effect  has 
not  been  noticed.  In  the  first  place,  we  should  recall  that  the  losses  will 
generally  be  rather  different  in  two  modes  even  when  they  have  the  same  frequency. 
As  a  result,  true  degeneracies  wall  be  relatively  rare  compared  with  the  milder 
quasi-degeneracies.  In  the  second  place,  we  see  that  even  the  quasi- 
degeneracies  are  not  particularly  common  in  the  small  thin  motors  typically 
used  in  research  studies  except  at  very  high  frequencies  where  one  would  expect 
the  general  damping  to  be  great  enough  to  keep  things  under  control. 

In  looking  for  this  phenomenon,  the  principal  clue  would  be  the  observation 
that  two  modes,  one  or  both  of  which  is  individually  stable,  together  become 
unstable  when  their  frequencies  cross. 

Price  has  observed  catastrophic  interaction  for  modes  having  the  same 
frequency  in  some  of  his  experiments.  It  appears,  however,  that  this  is  an 
interaction  between  certain  axial  and  tangential  modes.  In  linear  theory, 
axial  and  tangential  modes  arc  orthogonal  be  '  case  of  Jieir  different  angular- 
distribution,  and  thus  should  not  be'  coupled.  At  finite  amplitude,  however, 
in  the  presence  of  mean  flow  and  erosion,  this  exact  orthogonality  will  no 
longer  exist.  It  is  possible  that  Price's  observations  represent  the  same 
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phenomenon  where,  however,  the  coupling  arises  from  the  non-linearities 
implicit  in  the  large  amplitude  case. 

The  mode  maps  obtained  exp  riment^lly  by  Angelas  provide  an  appro¬ 
priate  dismay  of  data  for  which  a  number  of  crossing  points  of  mode  frequencies 
exist.  To  the  best  of  our  knowledge,  none  of  his  exp*  rimc~fs  have  shown  any 
particular  oscillations  attributable  to  these  intersections,  nor  have  they  shown 
the  catastrophic  behavior  observed  by  Price,  it  woe 'd  appear,  therefore,  either 
that  the  acoustic  erosion  constant  was  quite  small  in  i  'se  *  ~eqi  ncy  regions,  or 
that  th^  losses  in  the  intersecting  modes  were  such  that  noth  ng  approaching  a 
true  degeneracy  existed.  At  mis  time  there  no  lnmrmation  available  to 
decide  which  rl  these  reasons  was  dominant.  Wh  'ther  or  not  phenomena 
associated  with  mode  degeneracies  have  actuall;  occurred  in  researches  of  other 
investigators  is  essentially  not  determinable  mless  the  data  are  processed  ir.  a 
mannej  similar  to  that  used  by  An^elus. 

5.  3  Selection  Rules 

Up  ic  this  point,  we  have  not  considered  the  possibility  that  some  degenera¬ 
cies  (or  quasi -degeneracies  will  fail  to  contribute  because  of  zeros  in  the 
numerators  of  the  cr,ficai  terms  in  Eq.  (7).  Because  of  the  potentially  Urge 
degenerate  effects,  however,  the  selection  rules  introduced  by  zeros  in  the 
numerator  o  f  Eq.  1.7}  should  be  mt  ..tiuned. 

Referring  to  Eq.  (5e)  and  of  .  we  _  e  that  the  b.  and  c  .  will  vanish 
for  l  *  m  either  odd  or  even,  depending  on  symmetry  of  the  function  for  which 
they  are  the  Fourier  coefficients  Thus,  for  example,  analysis  snows  that  if 
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^  v  |  car;  oe  made  into  an  antisymmetric  function  about  the  midpoint  of  the 

chamber  by  adding  a  constant,  then  the  b„  vanish  for  £  +  m  even. 

f.  m 

Analogously,  if  |  can  be  made  symmetric  about  L/2  by  adding  a  constant, 
then  the  b ,  vanish  for  £  +  m  odd.  With  respect  to  the  c  .  ,  we  find  that  if 

Signum  (v  )  is  symmetric  about  z  =  L  2,  then  they  vanish  for  £  +  m  even  (or 
for  i  *  m  odd  if  antisymmetric),  so  that  we  should  find,  at  least  to  first  order, 
no  erosive  contributions  corresponding  to  such  "forbidden"  intersections. 

As  a  result,  the  selection  rules  arc  quite  different  in  the  end  nozzle  and 
center  nozzle  cases.  For  center  nozzle  rockets,  we  find  the  the  b’s  and  c's 
vanish  if  the  sum  of  their  indices  is  odd,  while  .or  the  end  nozzle  'mse,  these 
coefficients  vanish  if  the  sum  of  their  indices  is  e^en.  Thus,  in  center  nozzle 
rockets,  large  erosive  contributions  from  degenerrfe  configurations  become 
possible  only  at  those  configurations  which  were  forbidden  in  the  end  nozzle 


rockets  and  vice  versa. 


VI.  DISCUSSION 


The  initial  objective  of  this  study  was  merely  to  see  whether  acoustic  tueory 
would  indeed  confirm  the  r -it her  widely  held  interpretation  of  experiments  to  the 
effect  that  the  influence  of  erosion  on  stability  is  usually  a  miner  one  exc.  pt  for 
the  general  tendency  of  high  steady  erosion  to  stabilize.  We  have  met  with 
success  insofar  as  this  objective  is  ci-rccrned.  Eut  we  have  also  found  that 
under  certain  circumstances  which  are  presently  rather  except. unai,  rather 
than  common,  erosion  should  be  very  important  in  determining  the  answer 
to  the  stability  tiuestion. 
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For  the  most  part,  the  effects  of  erosive  velocity  have  been  separated 

into  two  components  which  arise  on  the  one  hand  from  the  erosive  component 

cf  the  me^n  flow  velocity,  and  on  the  other  hand  from  the  erosive  component 

of  the  acoustic  velocity.  With  respect  to  the  mean  flow  velocity  effects,  the 

analysis  indicates  that  if  the  dc  erosion  constant,  k,  and  the  mean  flow 

velocity,  v  ,  are  sufficiently  great  that  kf  v  |>a/1/3,  then  in  end  nozzle 
e  i  e  j 

rockets,  stability  should  be  significantly  enhanced  for  erosive  propellants  hav¬ 
ing  k>0,  and  significantly  diminished  for  propellants  having  k<0.  Experi¬ 
mental  evidence  does  substantiate  the  first  of  these  two  results,  but  we  are 
unaware  of  experimental  data  relating  to  negative  erosion  constants.  In  any 
case,  data  necessary  for  a  quantitative  comparison  with  experiment  do  not  yet 
exist.  The  odd  axial  modes  of  side  nozzle  rockets  should  be  somewhat  less 
sensitive  to  kv  in  linear  analysis.  We  know  of  no  quantitative  data  to  confirm 
or  deny  this  result,  except  again  to  note  the  fact  that  non- mesa  propellants 
in  side  nozzle  motors  seem  to  prefer  to  initiate  oscillation  in  odd  axial  modes. 

The  situation  with  respect  to  the  effect  of  the  erosive  component  of  the 
acoustic  velocity  is  numerically  somewhat  less  satisfactory.  The  analysis 
indicates  that  the  acoustic  erosion  constant  should  indeed  usually  have  little 
tearing  on  the  stability  question  -  except  in  two  kinds  of  exceptional  circum¬ 
stances.  Of  these  two,  the  first  pertains  to  the  side  nozzle  rocket,  and  the 
second  to  geometrical  configurations  for  which  the  system  becomes  acousti 
rally  degenerate. 

With  respect  to  the  side  nozzle  case,  we  have  found  that  acoustic  erosion 
contributes  in  first  order,  whereas  in  usual  end  nozzle  rockets,  this  effect 
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makes  its  first  appearance  in  the  second  order  terms.  Whether  acoustic 
erosion  increases  or  decreases  stability  (insofar  as  the  first  order  term  is 
concerned)  depends  on  *hs  imaginary  part  of  the  ac  erosion  constant,  K.  Thus, 
its  effect  must  vanish  at  sufficiently  low  frequencies  where  Kr**  k,  and  at 
sufficiently  high  frequencies  where  the  system  becomes  unable  to  respond  to  a 
sufficiently  rapidly  varying  erosive  velocity.  Unfortunately,  it  appears  that 
nothing  is  yet  known  regarding  the  relevant  frequency  interval. 

With  respect  to  the  degenerate  configurations,  a  realistic  quantitative 
calculation  in  the  immediate  neighborhood  of  these  ’’accidental"  degeneracies 
not  only  would  require  a  time -dependent  erosive  burning  theory,  but  also  a 
difficult  and  detailed  treatment  of  the  visco -elastic  losses  for  the  degenerate 
modes.  We  have  attempted  neither  of  these  tasks.  However,  an  order  of 
magnitude  calculation  has  shown  that,  at  least  in  principle,  the  effect  of  acoustic 
erosion  could  be  overwhelming  in  the  case  of  certain  degeneracies.  We  should 
perhaps  mention  a  third  atypical  configuration,  namely  rockets  having 
unusually  high  port  Mach  numbers.  Here,  too,  the  acoustic  erosive  velocity 
may  substantially  affect  the  balance  of  gains  and  losses,  with  the  sign  of  the 
effect  depending  on  the  phase  shift  associated  with  the  acoustic  erosion. 

Finally,  the  above  generalizations  must  be  tempered  by  the  thought  that 
for  some  future  generation  of  rockets  with  high  burning  rates,  highly  erosive 
propellants  and  short  fat  motors  many  features  not  previously  Observed 
may  well  become  apparent. 

With  respect  to  the  elementary  concepts,  and  to  accounting  for  the 
existence  of  the  phenomena  through  a  physical  picture,  one  may  now  easily 


# 
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visualize  their  occurrence  with  the  aid  of  a  little  hindsight.  The  effect  of  the 
dc  <  '•osion  constant  is  simplest  to  account  for,  so  we  shall  discuss  it  first. 

In  order  to  do  so,  then,  we  may  suppress  the  acoustic  erosion  by  setting 
K(cu)  =0  for  u)  fO.  Then,  referring  to  the  boundary  condition  at  the  burning 
surface,  we  see  that  the  turning  on  of  a  steady  uniform  erosive  velocity  merely 
increases  the  usual  steady  mass  flow  by  (1  +  k  |ve  |  ).  How  much  will  we  then 
have  to  alter  p/e  in  order  to  restore  neutral  stability?  Since  \i  was 
defined  as  the  fractional  increment  in  mass  flow  without  erosion,  we  would 
expect  to  find  the  increment  in 

just  as  is  given  by  Eq.  (6).  The  coefficient  b^  ^  "which  is  equal  to  unity  only 
for  the  constant  flow  case  just  considered,  may  be  regarded  as  describing  an 
appropriate  average  of  the  erosive  velocity  over  the  mode. 

The  possibility  of  a  large  acoustic  interaction  for  the  degenerate  cases 
might  also  have  been  anticipated.  For  simplicity,  let  us  now  suppress  the  dc 
erosion  by  setting  k  =  0.  For  an  end  nozzle  rocket,  any  single  mode  with  even 
l  has  the  acoustic  pressure  and  radial  velocity  symmetric  about  the  center  of 
the  chamber,  whereas  the  acoustic  axial  (erosive)  velocity  is  then  antisymmetric 
about  the  midpoint.  The  reverse  situation  Applies  if  1  is  odd.  Thus,  the  mass 
flow  perturbation  introduced  by  the  axial  velocity  tends  to  destroy  the  pressure 
distribution  which  gives  rise  to  it.  But  if  we  have  two  modes,  (at  the  same 
frequency),  one  with  even  and  one  with  odd  I,  then  the  erosive  perturbation  of 
one  mode  can  tend  to  reinforce  (or  cancel)  the  pressure  and  radial  velocity 
oscillations  of  the  other,  and  we  might  expect  a  relatively  large  effect.  For 
the  center  nozzle  rocket,  however t  the  antisymmetry  of  the  steady  state  flow 
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field  plays  a  dominant  role.  For  odd  modes,  the  pressure  distribution  is  anti 
symmetric  and  the  steady  flow  field  will  keep  the  velocity  distribution  sub¬ 
stantially  antisymmetric  also.  Thus,  large  interaction  between  the  erosion 
and  pressure  contributions  to  stability  would  be  expected  even  in  the  absence 
of  degeneracy.  On  the  other  hand,  if  the  mode  is  even,  the  situation  is  not 
essentially  different  from  the  end  nozzle  case.  Of  course,  the  dc  and  ac 
effects  are  not  always  so  neatly  separable,  so  that  this  backward  look  has 
qualitative  significance,  only.  But  it  should  probably  serve  to  caution  one 
against  any  ill-conceived  generalizations  of  these  results  to  other  configura¬ 
tions  with  new  boundary  conditions. 

In  conclusion,  it  should  probably  be  noted  that  the  present,  analysis  is 
hardly  a  ’’theory”,  but  rather,  a  mathematical  machine  designed  to  translate 
an  input  consisting  of  the  effect  of  acoustic  erosion  on  burning  rate  into  an  out 
put  consisting  of  the  resulting  effect  on  acoustic  stability.  We  have  investi¬ 
gated  the  output  of  this  machine  for  some  ad  hoc  inputs.  Further  progress 
requires  knowledge  of  the  inputs  which  may  be  encountered  in  reality.  This 
implies  either  a  time -dependent  theory  of  erosion,  or  direct  experimental 
measurement  of  the  dependence  of  burning  rate  on  fluctuating  gas  velocity. 

A  number  of  simplifying  assumptions  have  been  made  which  would  need 
examination  in  particular  cases.  The  cross  sectionally  uniform  flow  field 
assumption  is  certainly  incorrect,  but,  apart  from  its  axial  symmetry,  the 
results  depend  rather  insensitively  on  (space)  averages  of  v  ,  so  that 
perhaps  the  model  is  not  seriously  in  er*or  in  this  regard.  The  neglect  of 
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the  doppler  shifts  arising  from  the  mean  flow  field  may  become  significant 
for  high  port  Mach  numbers  (comparable  to  unity),  however,  and  the  present 
analysis  should  not  be  applied  in  such  instances.  Applicability  of  the  analysis 
is  further  limited  to  the  axial  boundary  conditions  which  we  have  treated  and 
it  is  not  too  clear  to  what  extent  the  results  would  be  modified  by  other  end 
boundaries.  Undoubtedly  the  selection  rules  would  be  quite  different,  in 
general*  The  neglect  of  acoustic  losses  in  the  gas  and  in  the  nozzle  is 
probably  most  significant  at  high  frequency  and,  in  the  case  of  end  nozzle 
rockets,  for  modes  with  high  axial  wave  numbers.  These  losses  are  signifi¬ 
cant  in  determining  which  frequencies  will  appear.  However,  we  have  been 
interested  in  the  relative  importance  of  the  pressure  and  velocity  effects,  so 
that  this  neglect  does  not  seem  to  be  a  significant  handicap. 

The  acoustic  boundary  condition  at  the  burning  surface  has  been  repre¬ 
sented  by  an  admittance  depending  on  frequency  but  not  otherwise  on  the 
detailed  nature  of  the  mode  being  considered.  Since  the  solid  itself  may  be 
a  participant  in  the  acoustic  motion,  this  representation  is  an  oversimplifica¬ 
tion,  The  solid  has  both  dilatation  and  shear  motion,  and  these  reflect  into 

a  component  (Y  )  of  the  surface  admittance  presented  to  the  gas  which  is  not 
s 

solely  frequency  dependent.  We  have  not  included  this  additional  realism  in 
the  calculations  presented  herein,  and  must  remind  the  reader  that  there  may 
be  occasion  upon  which  this  feature  may  be  quite  important. 
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